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Abstract 

Isotropic positive definite functions on spheres play important roles in spatial statistics, where they 
occur as the correlation functions of homogeneous random fields and star-shaped random parti- 
cles. In approximation theory, strictly positive definite functions serve as radial basis functions for 
interpolating scattered data on spherical domains. We review characterizations of positive defi- 
nite functions on spheres in terms of Gegenbauer expansions and apply them to dimension walks, 
where monotonicity properties of the Gegenbauer coefficients guarantee positive definiteness in 
higher dimensions. Subject to a natural support condition, isotropic positive definite functions on 
the Euclidean space M?, such as Askey's and Wendland's functions, allow for the direct substi- 
tution of the Euclidean distance by the great circle distance on a one-, two- or three-dimensional 
sphere, as opposed to the traditional approach, where the distances are transformed into each other. 
Completely monotone functions are positive definite on spheres of any dimension and provide rich 
parametric classes of such functions, including members of the powered exponential, Matem and 
generalized Cauchy families. The sine power family permits a continuous parameterization of 
the roughness of the sample paths of Gaussian processes on spheres. The paper closes with a set 
of sixteen research problems that provide challenges for future work in mathematical analysis, 
probability theory and spatial statistics. 
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1 Introduction 



For an integer A; > 2 let S'^"^ = {x G M'^ : ||x|| = 1} denote the unit sphere in the Euclidean space 
M'^, where we write for the Euclidean norm of x G M'^. A function : S'^ x S"^ — )• M is positive 
definite if 

n n 

^^QCj/i(xi,Xj) > (1) 

i=l 3=1 

for all finite systems of pairwise distinct points xi, . . . , x„ G S'^ and constants ci, . . . , c„ G M. 
A positive definite function h is strictly positive definite if the inequality in ([T]) is strict, unless 
Ci = ■ ■ ■ = c„ = 0, and is non-strictly positive definite otherwise. The function /i : S'^ x S'^ — )■ M 
is isotropic if there exists a function : [0, vr] — i- M such that 

h{x, y) = ^{9{x, y)) for all x,y e E>'^, (2) 

where 9{x, y) = arccos((x, y)) is the great circle or geodesic distance on S'^ and (■, ■) denotes the 
scalar product in M'^+^. Thus, an isotropic function depends on its arguments via the great circle 
distance 9{x, y) or, equivalently, the scalar product (x, y) only. 

For (i = 1, 2, ... we write '^d for the class of the continuous functions if) : [0, vr] — )• M with 
V'(O) = 1 which are such that the associated isotropic function h'm^h positive definite. Hence, 
we may identify the class with the correlation fu nctions of the mean-square continuous, sta- 
tionary and isotropic random fields on the sphere S'^ ( Jonesl [l963). We write for the class of 



the continuous functions ^ : [0, tt] — )• M with ^(0) = 1 such that the isotropic function h'm^h 
strictly positive definite, and we define = '^d \ ^1 ■ 

Recently, there has been renewed interest in the study of positive definite functions on spheres, 
motivated in part by applications in spatial statistics, where the members of t he class play 
crucial roles as the correlation functions of isotropic random fields on spheres (Baneriee', ^2005 



Huang et al. . 201l b. including the case of star-shaped random particles ( Hansen et al. , 2011, ). In 



a related development in approximation theory, the members of the class arise as radial basis 



functions for interpolating scattered data on spher ical domains (|Fasshauer and Schumakeii 11998 : 



Cavoretto and De RossiLl2010l : lLe Gia et al.Ll2010h . 



In this paper we study the classes \E'rf and with particular attention to the practically most 
relevant case of correlation functions and radial basis functions on a two-dimensional sphere, such 
as planet Earth. Applications abound, with atmospheric data assimilation and the reconstruction of 
the global temperature and green house gas concentration records being key examples. As it turns 
out, odd dimensions are mathematically more tractable, and so it is convenient to consider posi- 
tive definite functions on three-dimensional spheres, along with their natural restrictions to one- 
or two-dimensional spheres. Section |2] summarizes related results for isotropic positive definite 
functions on Euclidean spaces and studies their restrictions to spheres. In Section [3] we review 
characterizations of strictly and non-strictly positive definite functions on spheres in terms of ex- 
pansions in ultraspherical or Gegenbauer polynomials, and apply them to dimension walks, where 
monotonicity properties of the Gegenbauer coefficients in any given dimension guarantee positive 
definiteness in higher dimensions. 

The core of the paper is Section |4l where we supply easily applicable conditions for mem- 
bership in the classes '^d and and use them to derive rich parametric families of correlation 
functions and radial basis functions on spheres. In particular, compactly supported, isotropic pos- 
itive definite functions on the Euclidean space remain positive definite with the great circle 
distance on a one-, two- or three-dimensional sphere as argument. Consequently, compactly sup- 
ported radial basis functions on M?, such as Askey's and Wendland's functions, translate directly 
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Table 1: Parametric families of correlation functions and radial basis functions on one-, two- and 
three-dimensional spheres in terms of the great circle distance, 9 E [0, vr]. Here c is a scale or 
support parameter, r is a shape parameter, and a or z/ is a smoothness parameter, respectively, and 
we write t+ = max(t, 0) for t G M. The parameter range indicated guarantees membership in the 
classes ^j^, where d = 1,2 and 3. For details see SectionlH 



Family 


Analytic expression 


Parameter range 


Powered exponential 




c> 0; a G (0, 1] 


IVTatprn 


't^^^) r(u) \c) ^^{cJ 


r > 0- 1/ G CO -] 


Generalized Cauchy 




c> 0; a G (0,1]; r > 


Multiquadric 


(l-T)2^/(l+r2-2rcos0)'= 


c> 0; T G (0,1) 


Sine power 


ij{e) = 1- (sinf)" 


a G (0,2) 


Spherical 


m = {i + k^){i-^Jl 


c> 


Askey 


m = (1 - f ); 


c> 0; T > 2 


C2-Wendland 




c G (0,7r]; r > 4 


C4-Wendland 


^(0 = (l + rf + ^g)(l-f); 


c G (0,7r]; r > 6 



into locally supported radial basis functions on the sphere. A recent paper by iBeatson et al.l (|201 1|) 
provides criteria of Polya type that guarantee membership in the classes and we supplement 
these results in the case of spheres of dimension d < 7. Completely monotone functions are strictly 
positive definite on spheres of any dimension, and include members of the powered exponential, 
Matem and generalized Cauchy families. Some of the key results in the cases of one-, two- and 
three-dimensi onal spheres are sum marized in Table [H which complements a similar table in the 
recent work of Huang et al. ( 201 lb . 

The paper ends in Section[5l where we discuss challenges for future work from both theoretical 
and applied perspectives. A collection of sixteen research problems in mathematical analysis, 
probability theory and statistics offers challenges of diverse difficulty and scope. 



2 Isotropic positive definite functions on Euclidean spaces 

In this expository section, we review basic results about isotropic positive definite functions on 
Euclidean spaces and study their restrictions to isotropic positive definite functions on spheres. 

Recall that a function : M'^ x M"^ — )• M is positive definite if the inequality ([T]) holds for all 
finite systems of pairwise distinct points si, . . . , a;„ G M'^ and constants ci, . . . , c„ G M. It is strictly 
positive definite if the inequality is strict unless ci = ■ ■ ■ = c„ = 0. The function /i : M'^ x M'^ — )■ M 
is radial, spherically symmetric or isotropic if there exists a function : [0, oo) — )• M such that 

h{x,y) = ^{\\x -y\\) for all x,yeW^. (3) 



3 



For an integer > 1 we denote by the class of the continuous functions ip : [0, oo) — )■ M 
with ip(0) = 1 such that the function /i in © is positive definite. Thus, we may identify the class 
$d with the characteristic functions of spherically symmetric probability distributions, or with the 
cor relation func t ions o f mean-square continuous, stationary and isotropic random fields on M'^. 
SchoenbergI (|l938|) showed that a function Lp : [0, oo) — )■ M belongs to the class if and only 



if it is of the form 



/ T{d/2)(- ] J(d-2)/2(rt)dF(r) for t > 0, 



fit) 



(4) 



where J denotes a Bessel function (|Digital Library of Mathematical FunctionsLl201 iL Section 10.2) 
and F is a probability measure on [0, oo). The classes $d are nonincreasing in d. 



$1 D $2 D ■ ■ ■ D $c 



with the inclusions being strict, and if > 2 a ny noncons tant member (p of the class corresponds 
to a st rictly positive definite function on R*^ (ISunl . ll993l . Theorem 3.8). As shown by ISchoenberg 
(|1938b . the class $00 consists of the functions (p of the form 



exp(-rV)dF(r) for t > 0, 



(5) 



[0,00) 



where F is a probability measure on [0, 00). 

From the point of view of isotropic random fields, the asymptotic decay of the function at 
the origin determines the sm oothness of the associated Gaussian sample path almost surely (AdleA 
2009l : lGneiting et al.l . l2oTTb . Specifically, if 



^(0) - ^(t) = o{n 



as 



(6) 



for some a G (0,2], which we refer to as the fractal index, the graph of a Gaussian sample 
path has fractal or Hausdorff dimension D = d + 1 — ^. Ifa; = 2, the sample path is smooth 
and differentiable and its Hausdorff dimension, D, equals its topological dimension, d. If a G 
(0, 2), the sample path is non-differentiable. Table |2] shows some well k nown p arametric families 
within the class $ryi, namely the powered exponential family (Yagl oml 1987f). the Ma tern class 



(|Guttorp and GneitingL |2006|) and the generalized Cauchy family (|Gneiting and Schlatheii l2004|) . 
along with the corresponding fractal ind ices. Table |3] shows c ompactly s upported fam ilies within 



Jl9 



the class $3 that hav e been discussed by iWendlandl (|1995^) and iGneitingI (| 1 999al 12002 ) . 

Yadrenkd (|l983[) pointed out that if is a member of the class $fe for some k > 2, then the 



function defined by 



V'(^) = <^(2sinf) for 6 e [0, 



IT 



(7) 



corresponds to the restriction of the isotropic function h{xi,Xj) = ^{\\xi — Xj\\) from x M'^ to 



ys Ly'^ ^, with the chordal or Euclidean distance expressed in terms of the great circle distance, 
on the sphere S*^^^. Thus, given any nonconstant member ip of the class the function tp 



k-i- 



defined by (|7]) belongs to the class ^ 

Va ri ous authors have a rgued in favor of this constru ction , including iFasshauer and Schumaker 
(Il998h . IGneitingI (Il999ah . iNarcowich and Ward (l2002h and iBanerjeel (|2005|), as it readily gener- 
ates parametric families of isotropic, strictly positive definite functions on spheres and retains the 
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Table 2: Parametric families of membe rs of the class where denotes the rnodifie d Bessel 
function of the second kind of order v (IDigital Library of Mathematical FunctionsL I2OIIL Section 
10.2), c is a scale parameter, a and v are smoothness parameters, and r is a shape parameter, 
respectively. 



Family 



Analytic expression 



Parameter range 



Fractal index 



Powered exponential 




= exp(- 


-m 


0; 


a G (0, 2] 


a G (0, 2] 


Matem 




2'-'— 1 

- r(u) 




c> 0; 




min(2z^, 2) G (0, 2] 


Generalized Cauchy 




= (1 + 




c> 0; 


a G (0,2]; r > 


a G (0, 2] 



Table 3: Parametric families of compactly supported members of the class $3 with support param- 
eter c and shape parameter r. 



Family Analytic expression Parameter range Fractal index 



Spherical 


= 


+ 


cJ + 







1 


Askey 


ifit) = 


(1 - i)i 




0; 


r > 2 


1 


C^-Wendland 


Vit) = 






c> 0; 


r > 4 


2 


C^-Wendland 


Vit) = 


(1+-^ + ^ 




0; 


r > 6 


2 



interpretation of scale, support, shape and smoothness parameters. In particular, the mapping (Q 
from (y9 G $A: to G ^fe-i preserves the fractal index, in the sense that 

7/'(0) - ^(^) = as OiQ (8) 

if Lp has fractal index a G (0, 2], as defined in Q. Nevertheless, the approach is of limited flexibil- 
ity. For example, if G $3 it is readily seen from dH) that the function ip defined in (|7]) does not 
admit values less than inf^^o t~^sint = —0.2127 . . . Furthermore, the mapping 6^2 sin| in the 
argument of Lp in ([7]), while being essentially linear for small 0, is counter to spherical geometry 
for larger values of the great circle distance, and thus may result in physically unrealistic distor- 
tions. In this light, we turn to characterizations and constructions of positive definite functions that 
operate directly on a sphere. 



3 Isotropic positive definite functions on spheres 

Let d > 1 be an integer. As defined in Section [U the classes and consist of the functions 
•0 : [0, tt] — 7- M with ^(0) = 1 which are such that the isotropic function h'm^h positive definite 
and strictly positive definite, respectively. Furthermore, we consider the class '^'^ = \E'd \ of 
the non-strictly positive definite functions, and we define 

00 00 00 

d=l d=l (1=1 



5 



Schoenberg ( 1942 ^ noted that the classes and are convex, closed under products and closed 



under limits, provided the limit function is continuous. The classes and ^+ are convex and 
closed under products, but not under limits. 

We proceed to review Schoenberg's (1942) classical representation of the members of the 
classes and $00 in terms of ultraspherical or Gegenbauer expansions, along with the associated 
characterizations of strict positive definiteness in a more recent strand of literature in approxima- 
tion theory. Towards this end, we require some classical results on orthogonal polynomials (Digital 
Library of Mathematical Functions, 201 1, Section 18.3). Given A > and an integer n > 0, the 
function C^(cos 9) is defined by the expansion 



1 



^r"C^(cos0) for 9e[0,n] 



(9) 



n=0 



(1 + r2 - 2rcos^)- 

where r G (— 1, 1), and is called the ultraspherical or Gegenbauer polynomial of d egree n. For 
referen ce later on, we note that C^(l) = T{n + 2A)/(n! r(2A)). If A = 0, we follow ISchoenberg 
(Il942h and set 

C°(cos^) = cos{n9) for 9 E [0,7i]. 
By equation (3.42) of lAskev and FitchI (Il969h . 



[n/2] 

E 

k=0 



{n-2k + u)r{k + X - u)T{n - k + X) 
k\T{n-k + u + l) 



COS I 



(10) 



whenever X > v > and 9 E [0,n], where the coefficients on the right-hand side are strictly 
positive. This classical result of Gegenbauer will be used repeatedly in the sequel. 

The following theorem summarizes Schoenberg 's (1942) char acter ization of the cla sses \E'rf 
and \E'oo, along with the corresponding results of , Chen et al.l (l2003|) and lMenegattd (|1994h for the 
classes \E'^ and respectively. 

Theorem 1 (Schoenberg; Chen, Menegatto and Sun). Let d > I be an integer, 
(a) The class consists of the functions of the form 



(d-l)/2 



(cos^) 



n=0 



Md-l)/2 



(1) 



for 9E[0,n], 



(11) 



with nonnegative coefficients hn4 such that Y1^=q ^n,d = I. If d > 2, the class \E'^ consists 
of the functions in "^d with the coefficients bn,d being strictly positive for infinitely many even 
and infinitely many odd integers n. 



(b) The class \E'oo consists of the functions of the form 

00 

m = Y.^n{cos9r for 9E[0,7r] 



(12) 



n=0 



with nonnegative coefficients bn such that X^^o ~ ^- class ^+ consists of the 
functions in "if 00 with the coefficients bn being strictly positive for infinitely many even and 
infinitely many odd integers n. 
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If d = 1, Schoenberg's representation (fTTT) reduces to the general form, 



^{6) = J2 bn,i cos{n0) for 6 e [0, n] 



(13) 



n=0 



of a member of the class \E'i, that is, a standardized, continuous positive definite function on the 
circle. By basic Fourier calculus, 

bo,i = - I i'io) (14) 

and 

hni = - [ cosine) ij{9)d0 (15) 

for all integers n > 1. IfV^ € \l/^ , then & „ i is strictly positive for infinitely many even and infinitely 
many odd integers n (IMenegattol 1 1 995L Theorem 2.2), but the converse is not necessarily true, and 
a concise characterization of the class '^f remains elusive. If d = 2, the representation (fTTI) yields 
the general form, 



^P{e) = J2:;fflPnicos9) for ^e[0,7r] 



n=0 



(16) 



of a mernber of the class \l/2 in terms of the Legendre polynomial P„ of integer order n > 0. As 
BinghamI (1 1973b noted, the classes considered here are convex, and so (fTTI) . (fT2l) . (fT3]) and (fT6l) can 
be interpreted as Choquet representations. 

The general form (fT2l) of the members of the class ^oo can be interpreted as a power series in 
the variable cos 9 with nonnegative coefficients. For example, a standard Taylor expansion shows 
that the function defined by 



\2c 



m 



'1 + — 2r cos ( 



for 9 e [0, 



(17) 



admits such a representation if c > and r G (0, 1); furthermore, the coefficients are strictly 
positive, whence ^ belongs to the class ^+ . When d = 2, the special cases in (fTTI) with the 
parameter c fixed at ^ and | have been known as the inverse mul tiquadric and the Poisson spline. 



respectively (IFreeden et al.Lll997l : ICavoretto and De Rossil . 120101) . In this light, we refer to (fTTI) as 
the multiquadric family. 



The sine power function of ISoubeyrand et al.l (|2008[) . 



1 - 



9 

sm - 
2 



for 9 e [0, vr] 



(18) 



also admits the representation (fT2l) if a G (0, 2]. Moreover, the coefficients are strictly positive if 
a G (0, 2), whence ^ belongs to the class \1'+ . The parameter corresponds to the fractal index a 
in the relationship ([8]) and thus parameterizes the roughness of the sample paths of the associated 
isotropic Gaussian processes on spheres. 

The classes are nonincreasing in d > 1, with the following result revealing details about 
their structure. The statement in part (a) might be surprising, in that a function that is non-strictly 
positive definite on the sphere S'^ cannot be strictly positive definite on a lower-dimensional sphere, 
including the case §^ of the circle. 
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Corollary 1. 

(a) If i' E n "^d' for positive integers d and d', then tp G Similarly, if ip E \E'^ fl ^d' for 
positive integers d and d', then ip G \&^/. In particular, 



where the union is disjoint. 

(b) The classes are nonincreasing in d, 

with the inclusions being strict. 

(c) The classes are nonincreasing in d, 

with the inclusions being strict. 



^ oo ^ ^ 005 



(19) 



d>l 



d 1 



Proof. In part (a) we use an argument of lNarcowichI (| 19951) . while our key tool in proving that the 
inclusions in parts (b) and (c) are strict is Gegenbauer's relationship (flOl) . 

(a) \i d > d', it is trivially true that ip G \E'^ implies ip G \1'^. If d < d', suppose, for a 
contradiction, that ip ^ ^d'- By part (a) of Theorem [T] applied in dimension d' > 2, ip is 
either an arbitrary even function plus an odd polynomia l in cos 9, or an arb itrary odd function 
plus an even polynomial in cos^^. By Theorem 2.2 of iMenegattol ( 1995 ) we conclude that 
ijj G for the desired contradiction to the assumption that '0 G Thus, ^p G The 
proof of the second claim is analogous, and the statement in (fT9l ) then is immediate. 

(b) The inclusion 'i'd+i — trivially true. To see that the inclusion is strict, note that by part 
(a) of Theorem[T]and the relationship (flOl) any tjj G '^d+i ^drnits a representation of the form 
(fTTI) with bn,d being strictly positive for all integers n > 0. However there are members of 
the class \1'^ with bn,d = for at least one integer n > (|Chen et al.Ll2003b . which thus do 



not belong to the class ^ 



(c) The inclusion ^d+i C v]/^ is immediate from part (a). To demonstrate that the inclusion is 
strict, we show that if n > 2 then the function ip^O) = Cn ^^^'^{cosO) belongs to but 
not to '^'d-^-l■ For a contradiction, suppose that tp G ^d+i- Then by the relationship (flOl) . ^ is 
of the form (fTTI) with at least two distinct coefficients 6„ ^ being strictly positive, contrary to 
the fact that ip is an extreme point of the convex class ^d. D 

The following classical result is a consequence of Schoenberg's representation (fTTI) and the 
orthogonality properties of the ultraspherical or Gegenbauer polynomials. 

Corollary 2 (Schoenberg). Let d > 2 be an integer. A continuous function : [0, tt] — i- M with 
ipif)) = 1 belongs to the class if and only if 



Jn,d 



2n + d-liV{^)f 
23-'^7r T{d-l) 



(20) 



for all integers n > 0. Furthermore, the coefficient bn^ in the Gegenbauer expansion ( [77]) of a 
function -p E "ifd equals the above value. 
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Interesting related results include Theorem 4.1 of Nar cowich and Ward (l2002|) . which ex 
presses the Legendre coefficient bn,d of a function ^ of the form (|7]), where k = d+1, i n terms of the 



Fourie r transform of the respective member (p of the class ^d+i, and Theorem 6.2 of iLe Gia et al 
(|2010() . which deduces asymptotic estimates. 

We now provide formulas that express the Gegenbauer coefficient hn d+ y. in terms of the coeffi - 
cients bn^ and bn+2,d- In special cases, closely related results were used by iBeatson et al.l (|201 

Corollary 3. Consider the coefficients bn,d in the Gegenbauer expansion 4771) of the members of 
the class ^d- 



(a) It is true that 
and 

for all integers n>l. 

(b) If d>2, then 

bn,d+2 

for all integers n > 0. 



^0,3 — ^0,1 ~ 2^^'^ 



bn,3 = {bn,l - bn+2,l) 



(21) 
(22) 



{n + d-l){n + d) 
d{2n + d-l) 



(n + l)(n + 2) 
(2n + + 3) 



bn+2,a 



(23) 



Proof. We take advantage of well known recurrence relations for trigonometric functions and 
Gegenbauer polynomials. 



(a) The identity in (|2TI) holds true because 2 (sin ( 



-'0,3 



■ 1 - cos(20) for e e [0, 7r], so that 
cos(2^))V^(^) de = bo,i-^b2,i. 



To establish (l22l) . we note that 



C^(cos^)(sine)2 = sin^sin((n + 1)0) = - {cos{ne) - cos{{n + 2)6)) 

for 9 E [0, tt]. In view of the formulas (l20l) and (fTSl) for the Gegenbauer coefficients 6„ 3 and 
the Fourier cosine coefficients 6„ 1, respectively, the claim follows easily if we integrate the 
above equality over 6 E [0, vr]. 



(b) By equation (18.9.8) of Digital Library of Mathematical Functions! (|201 1|) . 



{n + d-l){n + d) 
(d-l) (2n + d + 1) 



+ l){n + 2) ^(,_i)/2 



(d-l) {2n + d + l) 



ri+2 



COS ( 



In view of the formula (|20l) for the Gegenbauer coefficients bn,d, we establish (1231) by inte- 
grating the above equality over E [0, tt] and consolidating coefficients. □ 
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If /c > 1 is an integer, so that d = 2k + 1 is odd, the recursion (|23]) allows us to express the 
Gegenbauer coefficient bn,d in terms of the Fourier cosine coefficients &„+2,i5 • • • , ^n+2fc,i- For 
instance, 

(n + 2)(n + 3) / n + 2 n + 1 \ 

&n,5 [b^,, - 2 ^ + ^ &„+4,l J (24) 

for all integers n > 0, where b^^ = 26o,i and i = if ^ > 1. Similarly, if d = 2k + 
2 is even, we can express the Gegenbauer coefficient bn,d in terms of the Legendre coefficients 
bn,2, bn+2,2, ■ ■ ■ , ^n+2fc,2- Furthermore, it is possible to relate the Legendre coefficients to the Fourier 
cosine coefficients, in that, subject to weak regularity conditions. 



k=0 



where 



Ci, 



(6*+2fc,l - ^n+2fc+2,l) 5 (25) 

1 • 3 • ■ ■ (2A; - 1) (n + 1) ■ (n + 2) ■ ■ ■ (n + fc) 
(2n)! fc! (2n + 3) • (2n + 5) • • • (2n + 2fc + 1) 

for integers n > and /c > 0, as used in special cases by Huang et"aD ( 2011 ). 

The recursions in Corollary [3] allow for dimension walks, where monotonicity properties of the 
Fourier cosine or Gegenbauer coefficients in a given dimension guarantee positive definiteness in 
higher dimensions, as follows. 

Corollary 4. Supose that the function ip : [O,^] ^ is continuous with ip{0) = 1. For integers 
d > 1 and n > 0, let bn,d denote the Fourier cosine and Gegenbauer coefficients 47?1) . 4771) and 
d20l) of ip, respectively. 

(a) The function ip belongs to the class ^3 if, and only if, 62,1 < 2 60,1 and 6n+2,i < bn^ifor all 
integers n > 1. It belongs to the class \E'^ if, and only if, furthermore, the inequality is strict 
for infinitely many even and infinitely many odd integers. 

(b) If d> 2, the function ip belongs to the class if and only if, 

2n + d + ?,{n + d-l){n + d) ^ 
- 2n + d-l in + l)in + 2) ^""^ ^^^^ 

for all integers n> 0. It belongs to the class \l'^+2 ^f' ^"^^^ tf' furthermore, the inequality 
is strict for infinitely many even and infinitely many odd integers. 

(c) If d>2, the function ip belongs to the class \1'^^_2 if bn+2,d < bn,dfor all integers n >0. 

Proof. Parts (a) and (b) are immediate from Corollaries [21 and [3l in concert with Theorem [H part 
(c) then follows from part (b). □ 



4 Criteria for positive definiteness and appfications 

In this section we provide easily applicable tests for membership in the classes and and 
apply them to construct rich parametric families of correlation functions and radial basis functions 
on spheres. 
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4.1 Locally supported strictly positive definite functions on two- and three- 
dimensional spheres 

There are huge computational savings in the interpolation of scattered data on spheres if the strictly 
positive definite function used as the correlation or radial basis function admits a simple closed 
form, and is supported on a spherical cap. Thus, various authors have sought such functions 
(|Schreineii 119971: iFasshauer and Schumakeii . 119981 : iNarcowich and Ward with particular 

emphasis on the practically most relevan t case of the two-dimensional sphere. 

The following result of iLevyl (1196 1|) constructs a positive definite function on the circle from 
a positive definite function on the real line, in the form of a compactly supported member of the 
class $1. Here and in the following, we write V5[o,7r] for the restric tion of a function Lp : [0, 00) — > M 
to t he inte rval [0, vr]. Essentially identical results can be found in lSasvaril (J 994 Exercise 1.10.25) 
andlWood (1995). 



Theorem 2 (Levy). Suppose that the function G $1 zi' such that ip{t) 
restriction ip = (Pio,tt] belongs to the class \E'i. 



for t > 71. Then the 



On the two-dimensional sphere, a particularly simple way of constructing a locally supported 
member of the corresponding class is to take a compactly supported member of the class $3, 
such as any of the functions in Table [3l and apply Yadrenko's recipe (|7]). If (/? G $3 is such 
that (f{t) = for t > a, where a G (0, tt), then the function ip in ^ belongs to the class 
and satisfies %jj{9) = for 9 > 2arcsin|. However, the approach is subject to the problems 
described in Section [2l in that the use of the chordal distance may result in physically unrealistic 
representations. 

Perhaps surprisingly, the following result shows that compactly supported, isotropic positive 
definite functions on the Euclidean space remain positive definite with the great circle distance 
on a one-, two- or three-dimensional sphere as argument. Thus, they can be used either with the 
chordal distance or with the great circle distance as argument. 

Theorems. Suppose that the functiomp G $3 is such that ip(t) = Ofort > vr. Then the restriction 
= '^[o,n] belongs to the class 

Proof. If the function G $3 is such that (p(t) = for t > tt, the Fourier cosine coefficient (fTSl) 
of its restriction ijj = ^[0,^] can be written as 

2 If 

bni = - / cos{n9)'ip{9) d9 = - cos{nt)^{t) dt = 2(j){n) 

for all integers n > 1. Here, the function : M — > M is defined by = (p{\t\) for t G M, 
and = 2^ /j^cos(Mt)0(t) dt denotes the inverse Fourier transform of evaluated at m G M. 
Similarly , 60 1 = 0(0). By equation (36) of GneitingI (1998a) in concert with the Paley -Wiener 
theorem (IPaley and Wienei . [1934), the function (j){u) is strictly decreasing in m > 0. Thus, the 
sequence bn,i is strictly decreasing in > 1; furthermore, 62,1 = 20(2) < 20(0) = 2 60,1 • By part 
(a) of Corollary m we conclude that G □ 

Theorem [3] permits us to use any of the functions in Table[3]with support parameter c < vr as a 
correlation function or radial basis function with the spherical or great circle distance on the sphere 
S'^ as argument, where d < 3. This gives rise to flexible parametric families of locally supported 
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members of the class where d < 3, that admit particularly simple closed form expressions. A 
first example is the spherical family, 

2 



m 



16 



1 - - 

C 



for ^G[0,7r], 



(27) 



which derive s from the family (l34l) in the c l ass that serves as a popular correlation model in 
geostatistics dOneiting , 1999ch . Huang et"aD ( 20 1 1 . Section 3.2) showed that the function in (ITT]) 
is in the class ^2 if c < tt; Theorem |3] yields the stronger result that tp E if c < n. 

Usi ng fam i lies of cornpactly supp orted members o f the class $3 introduced and discussed by 



Askeyl (Il973|) . IWendlandl (| 19951) and lOneitingl (|l999ah . Theorem |3] confirms that the truncated 
power function 

^jj{e) = (1-- ) for ee [0,7r] (28) 



belongs to the class if r > 2 and c G (0, tt), as shown previously by iBeatson et al.l (|201l() . If 
smoother functions are desired, the Wendland function 



m 



1 



e 

r- 

c 



e 

1 - - 

c 



for 9 G [0, tt] 



is in \1'3 if r > 4 and c G (0, tt); similarly, the function defined by 

-2 1 / 



m 



9 

l + T- 

C 



T 



for ^G[0,7r] 



(29) 



(30) 



belongs to the class if r > 6 and c G (0, vr). 

In atmospheric data assimilation, locally supported isotropic positive definite functions play 
important roles in the distance-dependent filtering of spatial covariance estimates on planet Earth 
(IHamill et al.L 1200 ll : iBuehner and Charronl . 120071) . The traditional construction in this literature 
relies on the Gaspari and Cohn (1999) function. 



{8t^ + 8t 



0<t< 



2' 



< t < 1, 



0, 



1 

2 

t > 1 



which belongs to the class $3, along with all functions of the form ip(t) = ifGcit/c), where c > 
is a constant. Yadrenko's construction ^ then yields the localization function 



^1(6*) = ^GC\ 



sm 



sm 



for ^G[0,7r], 



(31) 



which is a member of the class with support [0, c]. Theorem [3] suggests a natural alternative, in 
that, for every c G [0, vr], the function defined by 



MO) = VGc(l] for ^G[0,7r] 



(32) 



is also a member of the class with support [0,c]. Clearly, ip2{9) > 'ipi{9) for 9 G (0,c), as 
illustrated in Figure [H where c = |. This suggests that ^2 might be a more effective localization 
function than the traditional choice, ^1, in operational data assimilation. Similar comments m ight 
apply to covariance ta pers in spatial statistics, as proposed and applied by iFurrer et al.l (|2006|) and 
Kaufman et al.l (|2008[) . among others. 
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THETA 



Figure 1: The lGaspari and Cohn localization function with Euclidean or chordal distance as 

argument (^i; equation (|3TI)). and with spherical or great circle distance as argument (^2; equation 
(1321)). respectively, where c = |. While both functions are members of the class \E'^ and have 
support [0, c], we note that ^'2(6') > tpiiO) for 9 G (0, c), which suggests that ^jj2 might be a more 
effective localization function. 



4.2 Criteria of Polya type 



Criteria of Polya type provide simple sufficient conditions for positive definiteness by i mposing 
convexity conditions on a candidate function and/or its higher derivatives. iGneitingldlOOlb reviews 
and develops Polya criteria on Euclidean space s; in what follows, we are concerned with analogues 
on spheres that complement the recent work of lBeatson et al.l (12011 ). 



The following th eorem summarizes criteria developed by lWoodI (ll995|) . lGneitingl (|1998b|) and 
Beatson et al.l (|2011|) in the case of the circle . 



Theorem 4. Suppose that the function t/^ : [0, vr] — ?► M continuous, nonincreasing and convex 
with -0(0) = 1 and ipiO) dO > 0. Then ip belongs to the class \E'i. It belongs to the class 
unless it is piecewise linear. 

Our next result is a Polya criterion that applies on spheres S'^ of dimension d < 3. Its proof 
relies on the subsequent lemmas, the first of which is well known. 

Theorem 5. Suppose that ip : [0, 00) -^^is a continuous function with y9(0) = 1, lim^^oo ^if) = 
and a continuous derivative such that —i^'{t^^'^) is convex for t > 0. Then the restriction 
= V^Io.tt] belongs to the class \E'3 . 

To give an example, we may conclude from Theorem [5] that the truncated power function (1281) 
with shape parameter r > 2 belongs to the class ^j]" for all values of the support parameter c > 0, 
including the case c > vr, in which it is supported globally. 

Lemma 1. Let d > 1 be an integer Suppose that C is a subset of a Euclidean space, and let P be 
a Borel probability measure on C. If ipc belongs to the class for every c E C, then the function 
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ip defined by 

i,{e)= / Vc(^)dP(c) for [0,7r] 



c 

belongs to the class ^d, too. If furthermore ipc belongs to for every c in a set of positive 
P-measure, then ip belongs to the class 

Lemma 2. For all c > 0, the function defined by 



MO) = [^+2c)y~c) ^ ^ ^^^^ 

belongs to the class 

Proof. By Theorem[3l the function tpc belongs to the class if c < tt. If c > vr, equations (fT4l) 
and ([T?]) yield the Fourier cosine coefficients 

bo,i = ^ (Sc^^ - 6nc^ + n') and h,i = ^{4 + 2c^ - n^). 
For k = 1, 2, . . ., we find that 

37r2 1 

02fc,l 



^3 fc2 

and 



3 / 2c^-7r^ 2c^-7r^-2 \ 
^'^ ~ 2^V(2't + l)^ ~^ (2A; + 1)4 J 



Thus, 62,1 < 260,1 and 62A,+2,i < 62^,1 for all integers k > 1. As regards the odd Fourier cosine 
coefficients, let A; > 1 and note that 



h h - 12 k + 1 



where 

g{k) = {2c^ - n^){16k^ + AAk^ + 88k'^ + A8k + 9) + 8(4A;^ + 8A; + 5) > 0, 

whence 62^+3,1 < ^2fc+i,i for all integers k > 0. Therefore, we conclude from part (a) of Corollary 
IHthat ipc belongs to the class for c > n. □ 



Proof of Theorem\5\ By Theorem 3.1 of lGneitingI (Il999c|) . the function ip admits a representation 
of the form 



where (pc is defined by 



cp{t) = / cpc{t) dF{c) for t > 0, 

i(0,oo) 



and F is a probability measure on (0, 00). Therefore, the restriction 'ip = ^[o^tt] is of the form 



tPie)= / MO)dF{c) for ee[0,n], 
'(0,00) 
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where V^c is defined by (1331) . By LemmaO the function tpc belongs to the class \1'^ for all c > 0, 
and so we conclude from Lemma [T] that ^ is in ^^j^. □ 



Our next theorem is a slight generalization of th e key result inLB eatson et al.l(|201 ll) . In analogy 
to the corresponding results on Euclidean spaces in lGneitingl(|1999fi) . the case n = 1 yields a more 
concise, but less general, criterion than Theorem[51 

Theorem 6. Let 7i < 3be a positive integer. Suppose that : [0, oo) — > M ?'5 a continuous function 
with (/)(0) = 1, limt^oo V^(^) = and a derivative y?*^") of order n such that (— is convex 
for t > 0. Then the restriction ip = fio^w] belongs to the class '^2n+i- 



In a recent tour de force, iBeatson et al.1 (|2011|) demonstrate the following remarkable result 
about the truncated power function, which is the key to proving Theorem [61 

Lemma 3 (Beatson, zu Castell and Xu). Let n < 3 be a positive integer, and suppose that c e 
(0, tt). Then the function defined by 



n+l 



for 9 e [0, tt] 



(35) 



belongs to the class \1' 



2n+l- 



The next lemma concerns the case c > tt, in which the truncated power function (1351) is sup- 
ported globally. 

Lemma 4. For all integers n > 1 and all real numbers c> tt, the function ipcn in d331) belongs to 
the class ^i,. 



Proof. By a convolution argument due to Fuglede (IHjorth et al.Lll998L p. 272), the function ^ip^^^if)) - 
1 — ^ belongs to the class ^oo- As the class \l'oo is convex and closed under products, we see that 
if n > is an integer and c > tt then 



1 



TT \ "+1 



for 9 e [0, tt] 



whence ipc,n is in the class \l/oo, too. By a straightforward direct calculation, the Fourier cosine 
coefficients of ^ |Jr.,^ an d ipc,2 are strictly positive for all c > tt. Therefore, by Theorem 1 of 



Xu and Cheney! (119921) and the fact that the class "^^^ is closed under products, the function ipcn is 



in the class for all integers n > 1 and all c > vr. Using part (a) of Corollary [H we conclude 
that V'c.n G for all integers n > 1 and all c > tt. □ 



Proof of Theorem^ By Theorem 3.1 of lGneitingI (|l999c|) . the function ip admits a representation 
of the form 



(^,,„(t) dF(c) for t>0, 



(0,oo) 



where the function (y9c,n : [0, oo) — t- M is defined by 



9 

1 - - 

c 



n+l 



Vc,n{t) =yl- -\ for t > 

and F is a probability measure on (0, oo). Therefore, the restriction ip = ^[o^tt] is of the form 



m 



[ V'c,„(^) dF(c) for 9e[0,TT], 
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where tpc^n is defined by (|35]) . By Lemmas [3] and |4] the function tjjc,n belongs to the class "^tn+i for 
all c > 0. Therefore, we conclude from Lemma[T]that ^ is in '^tn+v too- 



□ 



iBeatson et al. ( 2011 ) conjectured that the statement of Lemma [3] holds for all integers n > 1. 
In view of Lemma IH if their conjecture is true, the statement of Theorem [6] holds for all integers 
n > 1, too, with the proof being unchanged. 



It is interesting to observe that Theorem[6]is a stronger result than Theorem 1 .3 of lBeatson et al. 



(|2011|) . in that it does not impose any support conditions on the candidate function ip. In contrast. 
Theorem 1.3 ofl Beatson et al. ( 2011 ) requires ^jj to be locally supported. Such an assumption is 
also made in the following result. 



Theorem 7. Let c G (0, tt), and let : [0,7[] ^ M. be a continuous function with tp{0) 
^[0) = Ofore>c. 



1 and 



(a) Let k > \ be an integer. Suppose that ip has derivatives of all orders up to 2k on (0, tt) and 
define 



If there exists an a > ^ such that 



du 



=6»2 



k-1 



is convex on (0, tt), then ip belongs to the class \l'i. 



(b) Let k > be an integer Suppose that ip has derivatives of all orders up to 2k + 1 on (0, vr) 
and define 



mm 



±) ,p(V-u) 



=e2 



If there exists an a>\ such that 



d9 J 



is convex on (0, tt), then ip belongs to the class 



Proof. We first prove part (a). Define the funct ion : [0 , oo) — »• M by Lp{t) = ip{t) for t E [0, tt] 
and (f{t) = for t > tt. By Proposition 2.1 of iGneitingI (1200 ih . where we set I = 0, (p admits a 
representation of the form 



[ ipk(-]dF{c) for t>0, 

J{0,oo) \Cj 



(36) 



where the function ip^ '■ [0, oo) -> M is a certain member of the class $i that is supported on the 
unit interval [0, 1], and F is a probability measure on (0, oo). However, as ipit) = for t > tt, 
the probability measure F must be concentrated on the interval (0, tt]. Therefore, the candidate 
function tp admits a representation of the form 



ipc{9)dF{c) for 9e[0,7t], 



(37) 



(O.tt] 
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where ipc denotes the restriction of the mapping t t—^ '^k{\) to the interval [0, vr]. By Theorem [2l 
the function ifjc belongs to the class for all c G (0, vr]. Thus, we conclude from Lemma[T]that if) 
belongs to the class \E'i. 



The proof of part (b) is identical, except that we set / = 1 in Proposition 2.1 of lGneitingI (1200 1|) . 
the function ipk in (l36l) is a certain member of the class <I>3 that is supported on the unit interval, 
and by Theorem[3]the function in (|37]) belongs to the class ^t, for all c G (0, tt]. Therefore, we 
conclude from Lemma [U that ^ belongs to the class . □ 

When /c = in part (b) of Theorem d the function z/^c in the representation (l37l) is of the 
specific form (|271) . and the resulting criterion reduces to a weaker version of Theorem [51 The 
stronger statement in Theorem[5]depends on the fact that G '^^ for all c > 0, including the case 
c > vr in which ijjc is supported globally. This property depends on the asymptotic linearity of the 
function ( pn in (l36l) at the origin, as reflected by its fractal index, a = 1. By Proposition 2.1 of 



Gneitingl(|2001|) . the function (pt in the crucial representation (|36l) has fractal index a = 2 if A; > 1, 
in both part (a) and part (b), and thus the condition that c G (0, vr] is essential for the membership 
of the corresponding function ifjc in any of the classes as will become evident in Section I4~4l 

4.3 Completely monotone functions 

A function : [0, oo) — )■ M is completely monotone if it possesses derivatives ip^'^'^ of all orders 
with {—lYLp^^\t) > for all integers A; > and all t > 0. Our next result shows that the 
restrictions of completely monotone functions belong to the class ^+ . 

Theorem 8. Suppose that the function (p : [0, oo) — > M is completely monotone with (f{0) = 1 and 
not constant. Then the restriction ip = <y5[o,7r] belongs to the class \1'+. 

Proof. Let a > and consider the truncated power function (1351) . By Lemma IH ipn/a,n belongs to 
the class \1'+ for all sufficiently large integers n. Hence, the function ipa defined by 



MO) = e-"' = lim i^n/aAO) for 9 G [0, vr] 



belongs to the class \&oo, too. By Theorem |5l ipa ^ ^3^, which in view of part (a) of Corollary [T] 
implies that tpa G ^J^. Now suppose that the function Lp : [0, 00) M satisfies the conditions of 
the theorem. Invoking Bernstein's Theorem, we see that (p admits a representation of the form 



ip{t) = / e"'** dF{a) for t > 0, 

J\0,oo) 



'[0,00) 

where F is a probability measure on [0, 00). Therefore, the restriction ip = v5[o,7r] is of the form 

ip{9)= I e-"^dF(a)= I ipa{0) dF{a) for 9e[0,7i]. 

J[0,oo) J[0,oo) 

Since cp is not constant, F has mass away from the origin. As tpa is in the class for all a > 0, 
we conclude from Lemma[I]that ijj belongs to the class ^+ . □ 



Miller and Samkd (|200lh present a wealth of examples of completely monotone functions, 



which thus can serve as correlation functions or radial basis functions on spheres of any dimension. 
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4.4 Necessary conditions 



The next result states conditions that prohibit membership in the class and therefore in any of 
the classes Heuristically, the common motif can be paraphrased as follows: 

If a positive definite function admits a certain degree of smoothness at the origin, it 
admits the same degree of smoothness everywhere. 

The conditions in parts (a) through (d) of the following theorem can be interpreted as formal de- 
scriptions of circumstances under which this overarching principle is violated. Part (a) corre sponds 
to a w ell-k nown result in the theory of characteristic functions, parts (b) and (c) are due t o I Wood 
(|l995|) and lGneitingI (|l998bf) . respectively, and part (d) rests on a result of lOevinatzl (|l959() . 

Theorem 9. Suppose that the function ip : [0, tt] — t- M can be represented as the restriction 
ip = 0[o,7r] of an even and continuous function : M — !■ M that satisfies any of the following 
conditions. 

(a) For some integer k > 1 the derivative (f)^'^^\0) exists, but (j) fails to be 2k times differentiable 
on (0, tt). 

(b) For some integer k > 1, the function (j) is 2k — 1 times continuously differentiable on [— tt, tt] 
with 0(2fe-i)(7r) ^ and |0(2fc-i)(t)| < far some h > and Pi > 0, and is 2k + 1 
times continuously differentiable on [— 7r,0) U (0,7r] with \(f)^'^''^^\t)\ < b2\t\^^'^ for some 
62 > and 132 E (1,2). 

(c) For some integer k > I, the function is 2k times differentiable on [— vr, vr] with (p^^^n) 7^ 
for some odd integer j, where 1 < j < 2k — 1. 

(d) The function is analytic and not of period 27i. 

Then ip does not belong to the class ^1. 

In particular, as we will see in the examples in the next section, if a function ip : [0, 00) — )■ M 
is strictly positive and strictly decreasing, admits a relationship of the form ^ for some a E (0, 2] 
at the origin, and furthermore is such that the restriction ijj = (p[o,w] belongs to the class then 
necessarily a < 1. 



4.5 Examples 



We now apply the criteria of Sections l43l and l4r4l to study parametric families of globally supported 
correlation functions and radial basis functions on spheres. In doing so, we think of a > or c > 
as a decay or support parameter, of « > or z/ > as a smoothness parameter, and of r > as a 
shape parameter. 

First we consider the families in Table [2l namely the powered exponential, Matem and gener- 
alized Cauchy classes. For the parameter values stated in Table [21 these families serve as isotropic 
correlation functions on Euclidean spaces, where the argument is the Euclidean distance. We now 
investigate whether they can serve as isotropic correlation functions with spherical or great circle 
distance on S'^ as a rgument, thereby supplementing and completing the corresponding results of 
Huang et al. (l201 lh in the case d = 2. 



In Examples [Hand [21 the results and proofs do not depend on the scale parameter c > 0, and so 
we only discuss the smoothness parameter a or v, respectively. 
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Example 1 (powered exponential family). The members of the powered exponential family are of 
the form 



il){0) = exp 



for 9 e [0, tt] 



(38) 



If a < 1, a straightforward application of Theorem [8] shows that ^ G \E'^. If a > 1, we see from 
parts (b) and (c) of Theorem [9] that tp does not belong to the class where part (b) applies when 
a e (1,2), and part (c) when a > 2, both using k = 1. Therefore, if a > 1 then ip does not belong 
to any of the classes ^rf. 

Example 2 (Matern family). The members of the Matern family can be written as 

V'W = FrT for ^^e[0,7r], (39) 



where i^j^ denotes the modified Bessel function of the second kind of order u, as defined in 



Digital Library of Mathematical FunctionsI ( 201 ll Section 10.2). If v = n + ^, where n > is an 
integer, then tp equals the product of exp(— 0/c) and a polynomial of degree n in 9. In particular, 
if z/ = i then ij{9) = exp(-6'/c) for 9 E [0, tt]. If z/ = |, then 



1 + Q exp (^-^^ for 9 E [0, vr] 



and if u 



I then 



m 



9 19' 

- + 3^)-Pl-- 



for ^G[0,7r]. 



Responding to questions in Sections 3.3.1 and 3.3.4 of lHuang et al.l(|2011|) . we now seek conditions 
under which the Matern function tp in (l39l) is strictly or non-strictly positive definite with spherical 
distance as argument. 

We first show t hat if u < ^ then tp bel ongs to the class Towards this end, we note 

from Theorem 5 of iMiller and Samkol J200lb that if < | the function (pi : R -)■ R defined by 
( ^1 (t) = exp(t)t'''K^{t) for t > is completely monotone. By Theorem 1 of iMiller and Samko 
hOOl S the function 0o : K ^ K defined by 0o(t) = exp{-t)(pi{t) = fK^it) for t > is 
completely monotone, too. In view of Theorem [8] this proves our claim. 

Next we demonstrate that if z/ > | then tp does not belong to the class ^i, and thus neither to 
any of t he classes "^d- If £ (^^^^^ this follows from an application of the estimates in Section 
10.31 of Digital Library of Mathematical FunctionsI (|201 lb to part (b) of Theorem[9l where k = 1, 
/?! = 2z/ — 1 and /32 = 2z/ — 3. If z/ > 1, we apply part (c) of Theorem|9l where k = 1. The details 
are tedious but straightforward, and we omit them. 

In our next example the results do not depend on the scale parameter c > and the shape 
parameter r > 0, and so we discuss the smoothness parameter a only. 

Example 3 (generalized Cauchy family). The members of the generalized Cauchy family are of 
the form 

/ / 9\°'\~^^°' 
tP{9) = 1 + - for 0G [0,7r]. (40) 



+ 



If a < 1, a straightforward application of Theorem [8] shows that tp is di member of the class \1/ 
However, if a; > 1 the function 'ip in (l40l) does not belong to the class \E'i. If a G (1,2), this is 
evident from part (b) of Theorem |9l where k = 1, 13i = a — 1 and = a — 3; if a > 2 we apply 
part (c) of Theorem [9l where also k = 1. 
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In the context of isotropic random fields and random particles, the smoothness properties of the 
associated random surface are governed by the behavior of the correlation function at the origin. 
Specifically, if a function E ^2 admits the relationship ([8]) for some a G (0, 2], the corresponding 
Gaussian sample path on the two-dimensional sphere has fractal or Hausdorf f dimension 3 — | 
almost surely, as shown and described more explicitly by Hansen et alJ (|201lb . In this sense, the 
results in the above examples are restrictive, in that the smoothness parameter needs to satisfy 
a < 1 or z/ < i, respectively. For parametric families of correlation functions that admit the full 



range of viable Hausdorff dimensions, we refer to the sine p ower familv (flSl) of ISoubeyrand et al. 



(l2008h and the convolution construction in Section 4.3 of iHansen et al.l (120111) . Alternatively, 



Yadrenko's construction (|7]) can be applied to any of the parametric families in Table [2l 

In our final example it is convenient to work with the shape parameter r > 1 and the inverse 
scale parameter r > 0, where small values correspond to slowly decaying correlations. 

Example 4 (Bessel family). The members of the Bessel family are analytic functions of the form 



2 \(^-2)/2 



for 6 e [0, tt] 



(41) 



where J denotes a Bessel function. Their relevance stems from Schoenberg's representation @ 
of the members of the class in terms of the function in (|4T]) . where r = d. As a consequence, 
the analytic continuation of ^ to a function : [0, 00) — )■ M belongs to the class $d if and only if 
T > d, irrespectively of the value of the scale parameter r > 0. 

We now investigate whether the function belongs to the class "^d, starting with the case 
r = 1, in which (|4TI) reduces to 

iP{9) = cos{r9) for 9e[0,n]. 

This function -ip can be extended to an analytic function </) : M — t- M with shortest period 27r/r. 
If r = 1, then ip belongs to the class \E';^; if r = n where n > 2 is an integer it belongs to the 
class \E']~ \ \E'2^, based on arguments in analogy to those in the proof of part (c) of Corollary [U For 
non-integer values of the scale parameter r > 0, part (d) of Theorem[9]implies that ^(6') = cos{r9) 
does not belong to the class \l'i. 

Similarly, we see from part (d) of Theorem |9] that if r > 1 then the function ip in (|4TI) does not 
belong to the class "^i, irrespectively of the value o f : r > 0. Thus, it does not belong to any of the 
classes \E'd, as shown in Section 3.3 of Huang et"aD ( 2011 ) in the particular case r = 3 and d = 2. 



5 Challenges for future work 

In this paper, we have reviewed and developed characterizations and constructions of positive 
definite functions on spheres, and we have applied them to provide rich parametric classes of such 
functions. 

Some of the key results in the practically most relevant cases of one-, two- and three-dimensional 
spheres are summarized in Table [B Whenever required, the closure properties of the classes of the 
positive definite or strictly positive definite functions offer additional flexibility. For example, 
while all entries in Table [T] yield nonnegative correlations only, we can easily model negative cor- 
relations, by using convex sums or products of an entry in the table with a Legendre function of 
the form il;{9) = Ci°'"^^^^(cos^)/Ci'^~^^^^(l), including the case ip{9) = cos9 that arises when 
n = 1, irrespectively of the dimension d. 
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An important construction pri nciple t hat we have only touched upon here is spherical con- 
volution, which has been used by IWoodI (|l995|) . ISchreinen (|l997[) . Estrade and IstasI (|2010|) and 
Hansen et al.1 (1201 1|) to construct positive definite functions on spheres. For example, if we con- 
sider an indicator function k{9) = 1(9 < r) on the two-dimensional sphere, where r G (0, |] is a 
cut-off parameter, spherical self-convolution yields the member of the class \l/2 given by 



■m 



( 



TT — arccos 



(cos6' 



CSC r 



cot^ r) 



7r (1 — cosr) 

— 2 cos r arccos (esc cos r CSC r — cot cot r)) < 2r) for ^^G[0,7r], 



as shown in the appendix of lTovchigrechko and Vakseii(|2001h . In the same setting. lEstrade and Istas 
develop a recursive formula that yields closed form expressions for the membe rs of the class 



tha t arise from the self-convolution of uniform kernels on c/-dimensional spheres. iHansen et al. 



(|201lb apply the convolution approach to derive a parametric family within the class that 
admits the full range of viable Hausdorff dimensions for the associated Gaussian and Gaussian- 
related random fields. 

Despite substantial advances in the study of positive definite and strictly positive definite func- 
tions on spheres, many interesting and important questions remain open. I will thus close the paper 
with a set of research problems, which aim to stimulate future research in mathematical analysis, 
probability theory and spatial statistics. The problems vary in scope and difficulty, range from 
harmonic analysis to statistical methodology, and include what appear to be tedious but routine 
questions, such as Problem [2l along with well known major challenges, such as Problem [T] and 
Problems [T4] through [T6l which have been under scrutiny for decades. 

The first eight problems are of an analytic character and concern the characterization and 
breadth of the classes and respectively. 

Problem 1. Provide a concise characterization of the class ^1/+ of the strictly positive definite 
functions on the circl e in terms o f their Fourier cosine coefficients. Despite recent progress by 
Chen et al.l (l2003b an d fsuni (EoOS), and related advances for strictly positive definite functions on 
the real line (IDerrienl . 12010). the challenge is substantial, and a persuasive solution remains elusive. 

Problem 2. For integers k>l and n > 0, find the coefficients a„,i, . . . , an,k in the general formula 



bn,2k+l — 0'n,kbn+2i,l 



for the Gegenbauer coefficient 6„.2fc+i of a function ip G '^2k+i in terms of its Fourier cosine 
coefficients h 



'n,l, On+2,1, 



, bn+2k,i, of which equations (|22l) and (1241) are the first two instances. 
Similarly, find the coefficients in the general formula for the Gegenbauer coefficient hn,2k+2 in 
terms of the Legendre coefficients 6„_2, bn+2,2, • • • , &n+2fc,2- 

We now let \1'^ denote the class of the functions ^ with il){6) = Q fox > c, where 
c? > 1 is an integer and c G (0, vr]. If c < vr, the corresponding correlation function or radial basis 
function has local support. Similarly, we let $^ denote the class of the functions (/? G $d with 
= for t>c, where (i > 1 is an integer and c > 0. 

Problem 3. Using the notation just introduced, \E'^ is the class of the functions if) E '^d with 
'(/'(vr) = 0, and $J is the class of the functions G $d such that (p{t) = if t > vr. Theorems [2] 
and [3] express the fact that 

$d c n (42) 
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when d = 1 and d = 3, respectively, and these inclusions are critical to the constructions in 
Section m Does the relationship (l42l) hold true for all integers d > 1? I have not been able to 
find a counterexample, nor have I made progress towards a proof. A positive answer would reap 
immediate benefits, such as far-reaching generalizations of Theorem |6] and Theorem [71 

Problem 4. For an integer d>l and c G (0, tt], find 



il = inf^gv^c [-tfj" 



(43) 



As hinted at in Section l4~n this is a problem of applied interest when d = 2. In atmospheric data 
assimilation, locally supported isotropic correlation functions are used for the distance-d ependent 



reduc t ion of global scale covarian ce estimates in ensemble Kalman filter settings ( Hamill et al. 



20011: iBuehner and Charronl . |2007|) . In this context, it is desirable to use a member of the class 
\l/2 with minimal curvature at the origin, as formalized by the optimization problem in (|43l) . By 
Theorem [2] and Theorem [31 



1 4 



(44) 



< inf^g^c [-v9"(0)] = j(d_2)/2 
1 and d = 3, respectively, where ji, denotes the first positive zero of the Bes sel function 



when d 

J„. T he right-hand equality in ([441) holds true for all integers d > 1 and all c > (|Ehm et al. 



2004 Section 5). 

Problem 5. For an integer d>l and c G (0, tt] 
c on the sphere S^. In other words, find 



find the Turdn constant of a spherical cap of radius 



where x is an arbitrary point on S'^ and the integral is with respect to the Haar me asure. When d = 
1, the probl e m has been resolved in a strand of liter ature that includes the works of IStechkinl (| 1 972|) . 
Gorbachev! (|2001|) and llvanov and Ivanovl (|2010l) . Extant general ization s to higher dimensions 
apply to convex bodies in Euclidean spaces or to the torus (|Reveszl . 1201 ih . but not to spheres, and 
so the case d > 2 remains open. 

Problem 6. For an integer > 1, c G (0, tt] and 9 G (0, c), find 



When d = 1, this pointwise version of Tur an's problem that has bee n studied by lArestov et al 



(|2003|) . iKolountzakis and Reveszl (|2006|) and llvanov and Ivanovl (|2010|) . I am not aware of any re- 
sults in dimension d> 2, except that lower bounds on b'^jO) become ava ilable from existing lower 
bounds in the analogous problem on Euclidean spaces (|Ehm et al.L |2004|) together with Theorem |3l 
and Yadrenko's construction 

Problem 1. What are the smoothness properties of the members of the class \E'rf? In the Euclidean 
case, the members of the class admit a derivative of order [{d — l)/2] on (0, oo), where \c\ 



denote the greatest integer less than or equal to c, and this is the best estimate possible (|GneitingL 
1999bl) . In the case of spheres, a natural conjecture is that the members of the class admit a 
derivative of order \{d — l)/2] on the open interval (0, tt), and that this is the best estimate possible. 
In this context, it is of interest to recall from the representation ([T2l) that the members of the class 
\E'oo admit an absolutely convergent power series expansion in the variable cos 9, and thus have 
derivatives of any order on (0, tt). 
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Problem 8. When does the isotropic positive definite function ^ that is repr esented by a member 
-0 of the class ad mit a spherical convolution representation in the sense of Schreinen (11997^) and 
Hansen et al.1 (|201l|) ? If a convolution root exists, can it be taken to be real-valued and isotropic? If 



for ^ > c, can we find a convolution root that is supported on a spherical cap of radius |? 



Ehm et alJ (|2004l) have answers to analogous questions for positive definite functions on Euclidean 
spaces, many of which are surprising. 

Next we state open questions about the parameter spaces for the truncated power family, the 
generalized sine power family, and the quartic exponential model used by the Berkeley Earth 
Surface Temperature project. 

Problem 9. It is tempting to conjecture that if c G (0, tt) then the truncated power function (1281) 
belongs to the class if an d only if r > + 1). \i d = 1, it is easy to see that the statement 
is true. iBeatson et al.l (120111) proved the sufficiency of the condition when d = 3,5, and 7, and 



conjectured it to be sufficient for all odd integers d > 1. Similarly, I conjecture that if c G (0, vr) 
then the locally supported functions (|29l) and (l30l) belong to the class if and only if r > ^(d+5) 
and T > ^{d + 9), respective l y. Note that if th e inclusion (l42l) holds tr ue in any given di mension d, 
then results of lAskevI ( fl973l) . lOolubovl (Il98lh . IWendland (1995.) and lGneitingr(ll999ah imply the 
sufficiency of the aforementioned conditions. 



Problem 10. A natural generalization of the sine power family (fTSl) of ISoubeyrand et al.l (|2008[) 
takes the form 

sin|-j 1(6 <c) for ee[0,7r]. (45) 

In addition to the smoothness parameter a G (0, 2], the generalized sine power family admits a 
scale or support parameter c > 0; if c < vr, it is locally supported. For what values of a G (0,2] 
and c > does the generalized sine power function (l45l) belong to the class ^'d? As the class 
does not admit locally supported members, the ans wer will depend deli catel y on the dimens ion 
d. The criteria in Section |4] apply, and the results in IPasenchenkd (| 19961) and iGneitingI (|2000h in 
concert with Yadrenko's construction (|7]) can provide sufficient conditions. 

Problem 1 1. The key component of the Berkeley Earth Surface Temperature project feohde et al. . 
201 ih correlation model on planet Earth is isotropic, depending on spherical distance as 



ip{9) = exp 




for 9 G [0, tt] 



For what values of the parameter vector (c^ , Co , c.cj, d )' does this function belong to the class $2? In 
particular, do the values fitted by iRohde et al.l (|201 iL page 11) provide a valid correlation function 
on the two-dimensional sphere? Unfortunately, the fitted function does not satisfy the assumptions 
of Theorem[51 nor does part (b) of Theorem |7] apply, as is globally supported. 

The next two problems relate to the fractal index and the sample path properties of Gaussian 
and Gaussian-related random fields. 



Problem 12. When d = 2, [Hansen et al.l (|201l|) prove and illustrate that if a function i/j G '^d has 
fractal index a G (0, 2], as reflected by the relationship ([8]), then any suitably defined, nontrivial 
Gaussian sample path on a rf-dimensional sphere has fractal or Hausdorff dimension c/ + 1 — |. 
Prove this result for spheres of any dimension d > 1. 

Problem 13. In the case of Euclidean spaces, Abelian and Tauberian theorer ns relate the behavior 
of a correlation function at the origin to that of its Fourier transform at infinity (|Leonenko and Olenko 
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199ll : ISteinL Il999f) . Find general analogues on the sphere that relate the behavior of a mem- 
ber of the class "ifd at the origin to the asymptotic decay of the Gegenbauer coefficients in 
(l20l) as n — 7- oo. In particular, what asymptotic behavior is associated with a given fractal index 
a G (0, 2]? In dimension d = 1, the broad principle is that the fractal index a G (0, 2) correspon ds 
to Fourier coefficients Ka decaying like n~°'~\, though the details are intricate (|BoasL ll967alJbl: 



to Jr'ouner coerncients 6r,,i decaymg iiKe r . 
Binmore and S tratto nl 1 19691: IWolfe. 'l973: 



results and tools of iMalyarenkoi (i2004|) and 



ChanLll985h. The recursio ns in Corollary |3] and the 



zu Castell and Filbid (l2005h seem relevant here. 



We close with challenges in spatial and spatio-temporal statistics, where non-isotropic and 
non-stationary correlation models on spheres are in strong demand, with the isotropic correlation 
functions discussed in this paper serving as building blocks for these more complex models. 

Problem 14. The class \E'd can be identified with the class of the correlation functions of the mean- 
square continuous, stationary and isotropic real-valued random fields on the sphere S''. However, in 
many cases of practical interest random fields are vector- valued, comprising, for example, multiple 
environmental, geophysical or meteorological or variables on planet Earth. It is thus of great 
interest to characterize the corresponding classes of matrix-valued isotropic correlation functions, 
find analogues of Cramer's theorem in the Euclidean case, and identify flexible parametric families 
of such functions. While Ya drenko's construction (|7]) generalizes readily and can be applied to the 
multivariate Matem family (|Gneiting et al.l bOlOh . or to other matrix-valued corr elation fun ctions 
on Euclidean spaces, approaches that operate directly on a sphere are preferable (|junL 1201 1|) . 

Problem 15. Processes observed o n spheres are frequen tly non-isotr opic, with the topography of 
planet Earth gnon et al.L |2006|) . tot al column ozone JStein', '20gi and error field s from globa 
numerical weather prediction models ( Weaver and Courti er. 200*3" ^AiietaL. 2002) being exam- 
ples. It is therefore desirable to develop non-isotropic stochastic process models on spherical 
domains, similar to the development of non- stationary and no n-isotropic correlation mode ls on 
Euclidean spaces, such as t he space deformatio n technique of iSampson and Guttord (Il992h. the 
stochastic volati l ity m odel of^ Huang et al. ( 2011 ) and parametric approaches ("Paci orek and Schervish . 



2006c ISchlatheii 12010 ). On spherical domains , several app roaches have emerged, including ram- 
ifications of t he space defo rmation technique (|Dasl 120001) . methods using expansions in spheri- 
cal harmonics ISteinI (|2007|), approaches based on differential operators and nested stochastic par- 
tial differential equations (|jun and SteinL l2008l : iBolin and LindgrenL 1201 ih. and the approach of 



KatzfussI (|201lL Section 4.3) that builds on the work of iPaciorek and SchervishI (120061) and covari- 



ance tapering. Typically observed deviations from isot ropy include dependenci es of covariance 
structures on latitudes, with theory being provided by Hitc zenko and SteinI (|2012h . 

Problem 16. Frequently, the temporal development of a process observed on a sphere is also of 
interest, so that the process needs to be modeled on the sphere cross time. Ne vertheless, the 
literature on the corresponding correlation structures is sparse, with the work of I Jun and Stein 
(E007n being a notable exception. 

The challenges in these problems are complementary, thus calling for the development of si- 
multaneously multivariate, non-stationary and spatio-temporal stochastic process models on the 
sphere, along with statistical methodology for model fitting and model choice. A promising 
recent dev elopr nent i n this direction is the stochastic partial differential equations approach of 



Lindgren et all (I20 111) . 
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